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ABSTRACT 

In this paper we consider a o? = 2 random Ising system on a square lattice with 
nearest neighbour interactions. The disorder is short range correlated and asymmetry 
between the vertical and the horizontal direction is admitted. More precisely, the vertical 
bonds are supposed to be non random while the horizontal bonds alternate: one row 
of all non random horizontal bonds is followed by one row where they are independent 
dichotomic random variables. We solve the model using an approximate approach that 
replace the quenched average with an annealed average under the constraint that the 
number of frustrated plaquettes is keep fixed and equals that of the true system. The 
surprising fact is that for some choices of the parameters of the model there are three 
second order phase transitions separating four different phases: antiferromagnetic, glassy- 
like, ferromagnetic and paramagnetic. 

PACS NUMBERS: 05.50.+q, 02.50.+S 
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1. Introduction 

Mean field spin glasses models have been studied and deeply understood both from a 
static and a dynamic point of view and key words like replica symmetry breaking, aging 
and ultrametricity, have become of very wide use in statistical mechanics of disordered 
systems The reason why spin glasses have attracted so much attention is probably 

more a consequence of the many successful applications to biological modeling (neural 
networks, immune system, adaptive evolution) then to their original scope limited to the 
description of disordered materials. For this reason and may be for objective technical 
difficulties most of the typical features which are very well established for the mean field 
models have not been found out for short range spin glasses. For example, it is commonly 
believed that a finite temperature glassy phase only exists for d > 3 spin glasses while in 
d = 2 one has only the paramagnetic phase. This is an almost surely true statement if 
one consider d = 2 spin system with independent bonds []6|-[TT[| and with vertical-horizontal 
symmetry but may be it is a false statement if one consider d = 2 spin asymmetric systems 
with correlated disorder. For example, in models with layered disorder the existence of 
a low temperature phase seems to be an established fact |1T^ - |1^ , nevertheless, one may 
think that these models are pathological since layered disorder is somehow a long range 
correlated disorder. 

In this paper we consider a d = 2 Ising system where there is both a short range 
correlation of the disorder and an asymmetry between vertical and horizontal direction. 
The specific interaction we chose is not motivated by a deep physical insight but it is 
merely dictated by technical reasons. Nevertheless, the model is not very artificial and 
the disorder correlation is limited to the fact that frustrated plaquettes always are present 
in near couples while the asymmetry only lies in a difference of strength of vertical and 
horizontal bonds. 

We solve the model using an approximate approach that replaces the quenched average 
with an annealed average under the constraint that the number of frustrated plaquettes is 
keep fixed and equals that of the true system. The surprising fact is that for some choices 
of the parameters of the models one can find four different phases. 

The paper is organized as follows. 

In section 2, after a brief overlook of the constrained annealing, we introduce our 
model with a particular attention to the concept of frustration; then we write the partition 
function with constrained frustration and the relative free energy. 
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In section 3 we derive the solution of the model, obtaining an expression for the free 
energy that can be computed via numerical methods. Moreover, the ground state energy 
is exactly found. 

In section 4 the conditions that yield to second order phase transitions are derived. 

In section 5 we describe the various behaviours of the model, showing a total of 
four distinct phases, three almost conventional (high temperature paramagnetic phase, 
ferromagnetic and antiferromagnetic phases at low temperature) and a fourth paramagnetic 
phase that we guess to have a 'glassy' nature. 

In section 6 we present our conclusions. 

2. Constrained annealing 

The model is defined on a square d = 2 lattice and the interaction is supposed to be 
effective only between nearest neighbours. The number of spins is = LM where M is 
the number of columns of the lattice and L is the number of rows. 

The vertical bonds are supposed to be non random and one can assume without loss of 
generality that they equal 1 while the horizontal bonds alternate; one row of all non random 
horizontal equal 1 bonds is followed by one row where they are independent dichotomic 
random variables which equal 1 with probability p and equal the negative value —7 with 
probability 1 — p (see fig. p. 

It follows that the Hamiltonian of our model can be written as: 

L M 

-f^iv = - ^ ^ {aijai+ij + Jijaijaij+i) (2.1) 
i=i j=i 

where aij = ±1 is the spin in the site located by the z-th row and the j-th column while 
the Jij are the horizontal bonds which equal 1 when i is even and are defined by 

{1 with probability 1 — p 

(2.2) 
— 7 with probability p 

when i is odd. 

The model is parameterized by 7 > and p and, in general, it is random, except in 
the two limit cases p ^ and p 1. In the first limit case p = all the couplings equals 
the unity and, therefore, we have the pure d = 2 Ising model ||T^. In the second limit case 
p = 1 the model is also not random, but while all the vertical couplings equal the unity, the 
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horizontal couplings alternate one row in which they are all positive and equal the unity 
to one row in which they are all negative and equal —7. In this second limit the model 
can be solved by standard transfer matrix methods [jl5| and it shows a low temperature 
magnetic phase; for 7 < 1 this low temperature phase is ferromagnetic while for 7 > 1 
there is horizontal antiferromagnetic order and vertical ferromagnetic order between the 
spins. For the sake of simplicity we will call hereafter this complicated magnetic phase 
simply the antiferromagnetic phase. Finally, the special choice p = 1, 7 = 1 corresponds 
the so-called 'fully frustrated model' [jl^which is also not random and has a transition 
only at T = 

In order to explain the nature of our approximation, let us first recall that the elemen- 
tary unit for frustration is the plaquette. If the product of the signs of the bonds around a 
plaquette is negative the plaquette is frustrated, otherwise, the plaquette is unfrustrated. 
In our Ising model, only the sign of the random variable Ji^j with odd i can be negative, 
therefore the two square plaquettes which share this bond are frustrated if this bond is 
effectively negative (e.g. 'a' plaquettes in fig. ^) and they are unfrustrated if it is positive 
('b' plaquettes). As a consequence of this definition of elementary frustration, we may 
define the total frustration of the system as the rate of frustrated plaquettes. In our 
model 

9 ^ ^ 1 _ r . 

i=ij=i ^ ' 

This quantity equals p in average, furthermore the strong law of large numbers assures 
that 4>N ^ P with probability 1 in the thermodynamic limit. 

We are far from being able to solve the quenched model, nevertheless we think that the 
qualitative behaviour of the system is captured by the above definition of total frustration 
( [jl^, for a more general definition see [|18|). Therefore, our proposal is to consider an 
annealed approximation where 4>n is constrained to coincide, in the thermodynamic limit, 
with the quenched total frustration p. This model corresponds to averaging Z only over 
the realizations of the disorder with total frustration p. We not only believe that the 
approximated model has the same qualitative features of the quenched one, but it is also in 
good quantitative agreement with it. In fact, our experience is that constrained annealing 
is a really powerful tool for estimating the free energy of disordered systems ||T9| - |2^ . We 
would like also to stress that the fixed frustration model can be also seen as an independent 
model where the bonds as well as the spins are allowed to arrange themselves in order to 
minimize the free energy provided they satisfy the global frustration constraint. 
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In order to obtain the free energy of the fixed frustration model we follow the general 
method ([0, We must first define the generalized partition function 



ZN{P,'y,^i) = J2e^p[~PHN + ^iN{(j)N - p)] (2.4) 
(J 

where /5 = ^ is the inverse temperature, and the average is over all realizations of the 
couplings Jij, than we obtain the free energy of the constrained annealed model as 

f{(3, 7) = - min lim In Z{f3, 7, (2.5) 

where the N (yo limit means that both M and L must tend to the same limit. In fact, 
the minimization over /x automatically selects the realizations of the disorder for which 
(pN = p in the thermodynamic limit. 



3. Solution of the model 



The generalized partition function is a sum of a product of randomly independent 
variables, therefore, we can write 



L M 



zn{(3, 7, /^) = XI n n 

cr 1=1 j = l 



1 - J, 



+ 7 



P) 



The average can be now easily performed, obtaining: 



ZN{f3,'y,n) = exp 



N 

-(ln(4p(l-p))-2(2p-l)/i + a; 



) 5Z^^P 









(3.1) 



where we have introduced the effective hamiltonian 

L M 

Hn 



i=l j=l 

The new effective horizontal bonds Ji are not random, they are all equal in the same row, 
and they alternate two possible values in different rows; in fact, one has Ji = 1 when i is 
even and Ji = when i is odd. The constants a and b are 



a = In 



cosh P— + 77 In 

V 2 2 p 



cosh ( /3 



,1 + 7 



1 , 1 



+ a In ■ 

2^2 p 



cosh (/3i±2-/. + ilni^) 

6 = In + 

cosh + ^ _ 1 In 



(3.2) 



It is possible to show that b is a monotonic decreasing function of jj, with — 27/3 < 6 < 2/5, 
so that we can directly use b as variational parameter in order to realize the minimum in 

The effective hamiltonian Hn is indeed associated to a pure 2d Ising model with 
unitary strength couplings along the vertical bonds, and with alternated rows of unitary 
and strength couplings. This model can be solved by trivially generalizing the Onsager 
solution and it is mapped into the problem of diagonalizing a collection of 2 x 2 matrices. 
In the thermodynamic limit N ^ oo the total free energy ( |2.5| ) reads 



/(A 7) 



mm 

b 



b 1-p 6 + 2/37 p 2(3-h 

4^(l-2p)-^ln(e -l)-^Me -1) + 



+ 



1 



dq\n X{q, b) 



(3.3) 



47r/? Jo 

where X{q,b) indicates the maximum eigenvalue in modulus of the product of the two 
matrices 

Tp{q) = exp [(3*{tz cosq + sing)] exp(-2/3r;j) 
Tb{q) = exp [/3*(r2 cos q + sin q)] exp{-br;,) 
where (3* = — lntanh/3, and r^, Tz are Pauli matrices. After some trivial algebra one gets: 



X{q,b) = t{q,b) + ./t{q,br-l 



(3.4) 



with: 



, 9 sinh6 „ , „^ sinh(2/3+6) 

t{q, b) = 2cos^q . - 2cosQCOsh2/3 — _ _ ^_ ^ + 



sinh2/? 



sinh^ 2/3 



+ 



cosh2 2/3 cosh(2/3+fe) + cosh(2/3-6) 



(3.5) 



sinh^ 2/3 

The minimum of (|3.3| ) is realized for b = b* and it is achieved by looking for the zero 
of its derivative. One has the self-consistent equation for b: 



1 + 



2p 



2(1 -P) 



2(3 -b , , -2/37-6 
e — 1 1 — e 



dt{q,b) 
db 



'tiq,br-l 







(3.6) 



b=b* 



When p — > the model has to reduce to the standard Ising model and in fact, the previous 
formula leads to b* —>■ 2/3, while in other limit case p —>■ 1 one has b* —2^(3. 
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Equation (|3.6| ) is an ordinary equation in b, nevertheless, it cannot be explicitly solved 
so that we are not able to give a compact expression for b* in terms of T, 7 and p. However, 
( p.6|) and then (|3.3|) can be numerically computed with the necessary precision in order 
to fully investigate the model. Furthermore, at T = while we don't have the complete 
solution of we are able to derive the leading terms of b* and to compute the ground 
state energy Uq. We find out that Uq has different linear behaviours in p depending on 7 

C/o = -2 + (1 - 

The T = entropy can only be computed numerically and it is shown in fig. ^ 
Unlike the energy Uq, Sq equals a constant function of p for any of the three choices of 7. 
It is always zero for 7 < 1 and 5*0 > for 7 = 1. In the case 7 > 1 one can show that the 
entropy becomes positive for p > p, where p < |, and So{p = |) — 0.01. 



4. Transitions 

Let us stress again that formulae ( p.3| ) - ( |3.6| ) represent the solution of the model, and, 
in principle, all the informations about it can be derived from them. Fortunately, even if 
we are unable to give an explicit expression of b* starting from ( p.6|) , we can easily obtain 
some analytic results. For instance, in this section we find the conditions that yield to a 
second order phase transition. 

Looking carefully at (|3.3| ) - ( |3.6|) , one can realize that the mechanism of the usual 
Onsager transition is preserved: the discontinuity occurs when b*, the zero of ( |3.6| ), nullifies 
also the argument of the square root in (|3.4|), i.e. when 



t(g,5*) = l (4.1) 

(the case t{q, 6*) = — 1 is not possible since t{q, b) > 1 for Wq and V6). In other terms, one 
has to find out the solution b* of a system of two equation, (|3.6|) and (|4.1|) . Obviously this 
solution can exist only for certain values of 7 and p. 



A direct inspection of (|3.5|) shows that ( |4.1j ) can be satisfied only for the specific 
choices g = or g = tt, and it determines the existence of two distinct transition lines in 
the p — T phase diagram at fixed 7, (see fig. ^ ). The first line ends on the p = (pure 
Ising model) axis at the Onsager critical temperature, so that in the following we will refer 
to this transition line as the ferromagnetic one. The second line exists only for 7 > 1 
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and it ends on the p = 1 axis in correspondence of the critical temperature separating 
the antiferromagnetic phase from the paramagnetic one; for this reason we will call it the 
antiferromagnetic line. 

After some trivial algebra (|4.1| ) reduces to 



sinh(2/? + b*) = ±2^^^^ (4.2) 
^ sinh2 2/3 

where the sign + correspond to the ferromagnetic line {q = 0), and the sign — to the 
antiferromagnetic one [q = n). 

In the limit case p = {b* = 2/3), the ( f4.2| ) recovers the well-known result sinh(2/3) = 1, 
while in the other limit case p = 1 [b* = —2(3'j) the transition is present when sinh(2/3(l — 
7)) = ~'^ sinW-^%3 ^ finite temperature when 7 > 1 and at zero temperature for the 

'fully frustrated model' (7 = 1). 

From a practical point of view, in order to compute numerically the transition lines 
which are showen in fig. ^, it is convenient to solve ( [4.1D with respect to 6* 



±arcsinh(2^^^) -2/3 (4.3) 
' sinh^ 2/3 ' 



Then, keeping 7 fixed and substituting b* into ( |3.6|) , one obtains the two transition lines 
p(T), which can be easily computed by standard numerical algorithms. 



5. Phases 

In the previous section we have seen that (^4.2| ) gives the known critical temperatures 
of the non random models {p = or p = 1). Another preliminary information about the 
behaviour of the model comes from the observation that the right hand side of ( [4.2| ) goes 
to zero in the limit T ^ 0, so that the ferromagnetic and the antiferromagnetic lines must 
coincide when they end on the T = axis. Studying the leading terms of ( p.6| ) close to 
T = 0, one finds that this coinciding point is at p = 1 for 7 = 1 and at p = | for 7 > 1. 

The full description of the different behaviours can be derived computing the transition 
lines in the p — T phase diagram at varying 7, as seen in the previous section. The following 
four scenarios listed in fig. |^, are obtained. 

For 7 < 1 (see fig. where 7 = 0.8) only the ferromagnetic line is present, separating 
two well-known phases: a ferromagnetic phase at low temperature, and a paramagnetic 
phase at high temperature, exactly as for the Onsager non random model. In fact the 
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antiferromagnetic random couplings are too weak with respect to the ferromagnetic ones, 
so that they are not able to change the structure of the phases of the pure model. 

When 7 = 1 (fig. ||b) the scenario is quite similar to the previous one, apart from the 
fact that the ferromagnetic line reaches the axis T = at p = 1, in correspondence of the 
T = transition of the fully frustrated model. Notice that the antiferromagnetic transition 
line is still absent. In fact, the antiferromagnetic couplings have the same strength of the 
others but, for any p < 1, their number is lower than the number of horizontal ferromagnetic 
couplings so that the ferromagnetic order prevails at low temperature. 

The most interesting situation corresponds to the choice 1 < 7 < 2 (fig. |]c, where 
7 = 1.2). First of all notice that both the transition lines are present. The first line starts 
on the p = axis and it ends at p = | on the T = axis and it delimitates the low 
temperature ferromagnetic region. The second line starts on the T = axis at p = | 
ending on the p = 1 axis, and it delimitates the antiferromagnetic phase. Outside this two 
regions there is a non-magnetic phase, but notice that this one has a narrow tongue dividing 
the magnetic regions and reaching the T = axis at p = |. As a consequence, if one fixes 
the probability between | < p < | + Sp{'y), where ^^(7) is a small but finite number 
depending on 7, one can observe three different second order phase transitions varying the 
temperature T. The transitions separate four phases; starting from low temperatures, the 
first is ferromagnetic, the third is antiferromagnetic, the fourth is an ordinary paramagnetic 
phase while the second is a low temperature paramagnetic phase. In fig. ^ it is shown the 
specific heat C as a function of T at fixed 7 = 1.2 and p = 0.82, computed starting from 
the numerical solution of ( p.3| ) and ( |3.6| ). In this case the specific heat C exhibits three 
distinct peaks next to T ~ 0.299, T ~ 0.349 and T ~ 0.545. Indeed we need to magnify 
the picture since it is necessary to compute ( |3.3D with a great precision in order to show a 
certain growth of C around its discontinuity. 

The appearance of a low temperature paramagnetic phase between the antiferro and 
the ferromagnetic ones represents an interesting peculiarity of this model. In particular, 
it happens at relatively low temperature and with an extremely narrow width. These are 
the main features that persuade us to guess a glassy nature for this paramagnetic phase. 
Moreover in our constrained annealed model, as seen at the end of section 3, the region 
at low temperature with an unphysical solution (negative zero temperature entropy Sq) 
do not reach the critical transition point (p = | , T = 0) where the 'glassy' paramagnetic 
phase ends. 
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The description of the different scenarios is completed with the case 7 > 2 (fig. 
where 7 = 8). The structure of the phase diagram is similar to the previous one with the 
difference that the narrow tongue between the ferro and the antiferro phases is suppressed. 
As a consequence for p < | we only have the ferro and the para phases while for p > | we 
only have the antiferro and the para phases. When 7 — > 00 the temperature of end point 
on the p = 1 axis goes to infinite. 

6. Conclusions 

The surprising feature of our model is that for some choices of the parameters 7 and p 
the magnetic phases are separated by a low temperature paramagnetic phase. We do not 
expect any long distance magnetic order in this phase i.e. <crijaij+k> = <fi,jCi+fc,j> = 
in the limit k ^ 00 but we expect that <aijaij^k>'^ > and <aijai^k,j>'^ > in the 
same limit. Our proposal is that these last two quantities properly characterize the glassy- 
like paramagnetic phase and they should vanish in the high temperature paramagnetic 
phase. Unfortunately, it is well known that in two dimensional models the computation of 
long distance correlation is a difficult task and, in our case, the computation also involves 
averages over the disorder making the situation even more complicated. We think that 
some work can be made in this direction but it will demand much technical effort so that 
our claim that the low temperature paramagnetic phase is somehow a glassy phase is, at 
this point, more a conjecture than an established fact. 

Two more questions remain to be answered. The first is the most relevant: what is 
the role of the annealed approximation in the qualitative features of the phase diagram? 
Or better, is the new phase a mere consequence of the annealed approximation? In this 
case, our fixed frustration model would have an interest in itself but it would not be a good 
approximation of the quenched one. We think that this question can only be answered by 
direct Montecarlo simulation. 

The second question is: in this model the frustrated plaquettes appears only in couples, 
what is the role of this special correlation? To be more specific, a model in which all vertical 
bonds equals the unity while the horizontal are independent random variables which take 
two possible values of opposite sign would have the same qualitative behaviour? We cannot 
approximate such a model by our fixed frustration technique so that also this question can 
only be answered by Montecarlo simulation, nevertheless, we are convinced that the answer 
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should be positive. In fact, the special nature of correlation between plaquettes is short 
ranged and there is not reason why it should affect the long range behavior of the system. 

In conclusion we would like to stress that in spite of the very partial results contained 
in this paper and of the many unsolved questions this work sheds some light on the very 
important point of the existence of a finite temperature glassy phase for d = 2 frustrated 
systems. In fact, in the most restrictive interpretation of our result we can still affirm that 
the fixed frustration annealed d = 2 model has low temperature glassy-like phase, while 
in the most generous one we can say that a true glassy phase exists for quenched random 
d = 2 systems. 
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Figure Captions 

Fig. 1. Zero temperature entropy Sq as function ofpat7<l,7 = l and 7 > 1. In the 
last case becomes positive for p > p, where p < |, and 5'o(p = |) ~ 0.01. 

Fig. 2. p — T phase diagram at different 7: a) 7 = 0.8; b) 7 = 1; c) 7 = 1.2; d) 7 = 3. 

Fig. 3. Specific heat C as function of T for 7 = 1.2 and p = 0.82. The first peak, located 
at T ~ 0.299, is referred to an antiferromagnetic transition, while the others, 
magnified in the boxes (T ~ 0.349 and T ~ 0.545, to ferromagnetic transitions. 
The circles in the boxes represent our numerical data. 

Fig. 4. A typical realization of the system. The full lines represent the +1 ferromag- 
netic bonds, while the dashed lines are the —7 antiferromagnetic bonds. The 'a' 
elementary plaquettes are frustrated, at difference with the 'b' plaquettes. 
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Fig.l 
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